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ABSTRACT 

We have measured the radial velocities of five 51 Peg-type stars and one star 
known to be constant in velocity. Our measurements, on 20 A centered at 3947 
A, were conventional using Th/Ar comparison spectra taken every 20 or 40 min- 
utes between the stellar exposures. Existing IRAF routines were used for the 
reduction. We find grv < 20 m s _1 , provided 4 measurements (out of 72) with 
residuals > 5a RV are neglected. The observations were made with the CFHT 
Gecko spectrograph (R~110,000), fiber-fed by the CAFE system, a R y <10 m 
s _1 seems possible with additional care. This study was incidental to the main 
observing program and is certainly not exhaustive but the small value of cr RV 
implies that the fiber feed/image sheer system on Gecko + CAFE, essentially 
eliminates the long standing problem of guiding errors in radial velocity mea- 
surements. We are not promoting this conventional approach for serious Doppler 
planet searches (especially with Gecko which has such a small multiplex gain), 
but the precision is valuable for observations made in spectral regions remote 
from telluric lines or captive-gas fiducials. Instrument builders might consider 
the advantages of the CAFE optics which incorporate agitation and invert the 
object and pupil to illuminate the slit and grating, respectively, in future spec- 
trograph designs. 

Subject headings: instrumentation: spectrographs; techniques: spectroscopic, 
radial velocities; stars: planetary systems 
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1. Introduction 

Radial velocity measurements are of fundamental importance in astronomy. For stars, 
absolute radial velocities can be uncertain by as much as a km s" 1 because of motions 
in the stellar photosphere. Differential velocities, on the other hand, can be much more 
precise yielding important information about pulsation for single stars and masses for those 
in double or multiple systems. The most spectacular return from high precision radial 
velocities in the last decade has been the detection of the minute reflex accelerations caused 
by unseen planetary companions for some 100 solar-type stars (see, for example, Marcy & 
Butler (2000)). 

The basic design of spectrographs most commonly used to determine stellar Doppler 
shifts has changed little in over a century. The telescope images a star on the spectrograph 
slit which is collimated, dispersed and then reimaged onto a detector. The principal 
improvements have been in detectors, from eye to photographic plate to photomultiplier or 
solid state array. In modern spectrographs gratings largely replace prisms for dispersion. 
Precision in estimating wavelength displacements is compromised by motion and defocus 
of the star image at the slit. Both translate into a displacement of the stellar spectrum 
relative to the comparison spectrum. Flexure of the spectrograph decollimates Cassegrain 
instruments while, for bench-mounted spectrographs, flexure of the telescope causes its 
optical axis to wander in the spectrograph. Both effects lead to spectral shifts. 

Petrie & Fletcher (1967) thoroughly studied these effects with the McKellar 
spectrograph at the coude focus of the DAO 1.2-m telescope from photographic 
spectrograms. Working at a resolution of ~55,000, they found that even for bright stars 
with well modulated spectra, there were external errors of nearly 300 m s _1 , while for 
the sky they were closer to 100 m s -1 , and internal errors were even lower at ~70 m s _1 . 
They pointed out that without some scrambling of starlight at the slit, such systematic 
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errors would remain. The subsequent introduction of pupil or image slicers reduced this 
discrepancy by improving the uniformity of slit illumination. Campbell & Walker (1979), 
also used the McKellar spectrograph to demonstrate that, with a captive-gas to impose 
wavelength fiducials directly in the stellar spectrum, systematic errors could be dramatically 
reduced. Nowadays, with iodine vapor cells, errors have been reduced to a few m s -1 , the 
level at which precision begins to be limited by the natural velocity noise of the star's 
atmosphere (Marcy & Butler 2000). 

For many programs it is not always possible or appropriate to include either suitable 
telluric lines (generally water vapor) or captive-gas lines such as I 2 . This is certainly true 
for the blue/ultra- violet at high spectral resolution. In this paper we briefly demonstrate 
that, with minimal precautions, a precision of <20 m s~ x is possible for R>100,000 with 
a fiber-fed spectrograph over runs lasting at least several nights and probably for much 
longer. 

We became aware of this high precision following an observing run at CFHT in 2001 to 
observe Ca II H & K in the spectra of 51 Peg- type stars. These stars had been discovered 
to have Jupiter-mass secondaries with orbital periods of a few days from precise radial 
velocity measurements, mostly using captive gas fiducials. Such high quality radial velocity 
standards were not available to Petrie & Fletcher (1967) but offered us an excellent test of 
the Gecko + CAFE (CAssegrain Fiber Environment) radial velocity precision. In a run 
one year later we deliberately took comparison arc spectra before and after every stellar 
spectrum and it is these spectra which we discuss here. 

Our data were certainly adequate to calculate improved periods for each of the stars. 
The revised periods and phases will be published in another paper. 
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2. The Observations 
2.1. The spectra 

The observations were made on 26 to 30 July 2002 UT with the Gecko echellette 
spectrograph fiber fed by CAFE (Baudrand & Vitry 2000) from the Cassegrain focus of the 
Canada France Hawaii 3.6-m telescope (CFHT). Spectra were centered at 3947 A in the 
14th order which was isolated by a UV grism with some 60 A intercepted by the CCD. The 
dispersion was 0.0136A pixel" 1 and the 2.64 pixel FWHM of the thorium-argon (Th/Ar) 
lines corresponded to R= 110,000. The spectra were part of a long term program to monitor 
variations in the Ca II H & K reversals of stars with short periods (3 to 4.5 days) planets. 
For this reason, they were of high signal-to-noise in the continuum. The detector was a 
back-illuminated EEV CCD (13.5 /im 2 pixels) with spectral dispersion along the rows of the 
device. Single Th/Ar arcs were taken immediately before and after each stellar spectrum 
and, like the flat fields, were fed through the same fiber as the starlight. Probably more 
comparison spectra should have been taken but, as time was of the essence for this program, 
only single comparison spectra were taken. 

Fiber modal noise was suppressed by continuously agitating the fiber close to its output 
(see Baudrand & Walker (2000)). A Fabry lens at the fiber output projects the pupil as 
input to the spectrograph while simultaneously illuminating the grating with the object 
field. Baudrand et al. (1998) have pointed out that inversion of object and pupil improves 
both spectral resolution and spectral stability. In this arrangement the image projected on 
the grating is the fiber core output aperture enlarged 1000 times, while the fiber far field 
image is transformed into a pseudo slit by a four slice Bowen-Wallraven sheer of silicon 
which projects to some 50 pixels length in the spectrum. For this program there was no 
on-chip binning of pixels. 
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Because the CCD has a number of pixels that suffer from non-linear dark signal, we 
took a large number of dark exposures with integration times matching those of the various 
stellar, Th/Ar, and flat-field exposures, and created an average dark for each exposure 
time. Then, rather than using conventional biases to remove the baseline from each 
observation in the data reduction, the appropriate mean darks were subtracted from the 
stellar, Th/Ar, and flat-field exposures. Flat-fields were then normalized to a mean value of 
unity along each row and the exposures for all stars observed on a night were combined into 
a mean object to define a single aperture for the extraction of all stellar and comparison 
exposures, including subtraction of residual background between spectral orders (prior to 
flat-fielding). This aperture was ultimately used to extract one dimensional spectra of the 
individual stellar and comparison exposures and a one dimensional extraction of the mean, 
normalized flat-field (without an interorder background subtraction). The extracted stellar 
and comparison exposures were then divided by the one dimensional flat-field to obtain 
flat-fielded spectra. 

A specimen, flat-fielded spectrum of r Ceti is shown in Figure 1 with the ~ 20 A of the 
spectrum used to measure the radial velocities indicated. 

2.2. The radial velocities 

Table 1 lists the five program stars and the standard, r Ceti, plus spectral types and 
U magnitudes, the planetary orbital periods and velocity amplitudes, with details of the 
observations such as average S/N and exposure time. 

Radial velocities were estimated with the fxcor routine in IRAF 1 . The Th/Ar 

X IRAF is distributed by the National Optical Astronomy Observatories, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc. (AURA) under 
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comparison lines were used to provide dispersion-corrected stellar spectra. A Fourier 
cross-correlation was carried out on the dispersion-corrected spectra taking the first 
spectrum in the series for each star as the template. Hence, all differential radial velocities 
(ARV) are relative to the first spectrum on the first night. Both the template and the 
spectrum being measured were normalized with a low order polynomial and the correlation 
was taken over that part of the spectrum bounded by (and including) two strong aluminum 
lines (-3942-3963 A). 

In Figure 2 the ARVs for the five 51 Peg stars are plotted as a function of relative 
phase, and for r Ceti as a function of time. The sine curves have the published planetary 
orbital period and K max from Table 1 and were shifted in phase and ARV to give the best 
fit. For r Ceti there is a best fit is to a line of constant velocity. Below each velocity curve 
are shown the ARV residuals from the curves. Note that, because the amplitude of the 
stellar reflex velocity is very different in each case, the individual ARV scales vary widely. 
Values of o RV corresponding to the residuals are listed in Table 1. If four out of the 72 data 
points with a^v >100 m s" 1 are omitted, all of the a^v are <20 m s _1 and the two stars 
observed at the highest S/N have a^y ~10 m s -1 . 

2.3. The four extra-large residuals 

The motions of the comparison spectra expressed as ARV on each of the five nights 
are shown in Figure 3. The results are surprising. The ARV scales on nights 3 and 5 are 
an order of magnitude greater than on the other three nights. Nights 1 and 4 show a trend 
and modest scatter, with RMS about the trend being 9 and 18 m s" 1 , respectively. Night 2 
(only partially clear) shows no trend and an RMS scatter of 9 m s _1 . Nights 3 and 5 have 
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Table 1. The Program Stars 



star 


spectrum 


U 


exp time 


S/N a 


n b 


K c 


&RV 


P orb 








S 






m S _1 


m s _1 


days 


v And 


F7 V 


4.69 


1200 


490 


8 


70 


11 


4.6170° 


t Cet 


G8 V 


4.43 


1200 


620 


10 




9 




r Boo 


F7 VI 


5.02 


1200 


470 


12 


467 


31 (18) f 


3.3128 


HD 179949 


F8 V 


6.83 


2400 


250 


13 


102 


18 


3.093 


51 Peg 


G2 IV 


6.36 


1200 


270 


17 


56 


26 (18) f 


4.2293 


HD 209458 


GO V 


8.38 


2400 


150 


12 


87 


32 (18) s 


3.524738 



a average per 0.0136A pixel in the continuum 
b number of spectra 
c published values 

published orbital period: r Boo and v And (Butler ct al. 1997), HD 179949 (Tinney 
et al. 2000), 51 Peg (Marcy et al. 1996), HD 209458 (Charbonneau et al. 1999) 
°shortest of three planetary periods 

f value in brackets omits a single outlying point - see text 
s value in brackets omits two outlying points - see text 
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RMS scatters of 96 and 83 m s 1 , respectively! 

The most likely reason for spectral shifts on the detector are liquid nitrogen boil-off 
from the CCD Dewar, distortion of the coude room floor (Earth tides, dome stress etc.), 
and spectrograph 'seeing' effects. If these changes are slow and linear as they appear 
to be on nights 1, 2 and 4, then interpolation of the comparison line positions should 
properly calibrate the shifts for the stellar spectra. Most of the shifts are remarkably small 
considering that ARV=10 m s -1 corresponds to 1.3xlCT 4 A at 3947 A which is 10~ 2 of a 
pixel or 0.13 /xm. Monitoring the exposure meter to estimate the centroid of the stellar 
exposures would have allowed more accurate barycentric velocity corrections but we did not 
have such a system to hand. 

One might expect that the large, erratic residuals on nights 3 and 5 would be reflected 
in large reflex residuals for the stellar radial velocities on those nights but this is not always 
so. The delinquent velocity for 51 Peg on night 3 (see Figure 2) corresponds to the very 
large comparison residual (>300 m s -1 ). By contrast, the velocity for v And observed 
immediately afterwards which uses the same comparison spectrum agrees within a few m 
s _1 with the second value for v And on that same night. 

On the face of it, there is no consistent reason to reject the four large residuals. They 
may represent 'creaking' in the system which introduces occasional uncalibratable jumps. 
Unfortunately, because of the nature of the program we were not able to take comparison 
spectra more often than every 20 or 40 minutes. In two cases, the comparison lines show a 
residual of more than 100 m s" 1 which is reflected as an almost equally large residual in the 
stellar velocity. In these cases (51 Peg and HD 209458), using a mean Th/Ar which ignores 
the bad points restores the stellar velocities to the curve. In the other two cases (r Boo and 
HD 209458) there is no problem with the Th/Ar comparison. 

Some of the large residuals may be inherent to the reductions and not the spectrograph. 
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We used the same version of IRAF on two different platforms, PC-IRAF V2.11.3 and 
SUN/IRAF V2.11.3. The former generated several very large residuals in the Th/Ar 
comparison positions which disappeared on reduction on a SUN. The data presented in this 
paper were all reduced on a SUN but it is still possible that some of the large residuals are 
artifacts of the reduction. In short, we cannot be sure that occasional large residuals are 
not a risk with such conventionally determined velocities without a more extensive study. 

Differential variations between the strenghts of comparison arc lines can also limit 
precision. The temperature of the arc, its age, 'on' time, voltage stability etc., cause relative 
line strength variations and introduce line shifts especially when lines are unresolved blends. 
To judge from the low scatter in the comparison spectral shifts and values of Urv ~ 10 
m s" 1 for the two stars observed at high S/N, variations in the arc lines has not been 
important above this level and this grv may be the best achievable for single spectra with 
this spectrograph, but a more careful study is necessary. 



3. Conclusions 

The effectiveness of fiber agitation in overcoming modal noise has already been 
demonstrated (Baudrand & Walker 2000). In this paper we believe we have demonstrated 
that the optics of the CAFE fiber system which invert the conventional object and pupil 
(Baudrand et al. 1998) illumination of slit and grating together with agitation effectively 
eliminates the large guiding and tracking errors which once plagued conventional radial 
velocity measurements. Our experience suggests that comparison arcs taken with sufficient 
frequency are adequate to track slow spectrograph drifts for a bench mounted fiber-fed 
instrument in a stable environment. The precision depends on the resolution which, at 
110,000 for Gecko, is better than 10 m s -1 . Being an echellette format, the spectral 
curvature is small, but an echelle format would provide a multiplex gain of one or two 
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orders of magnitude over the single spectra captured by Gecko (in this paper we only used 
20 A!) . From the point of view of radial velocities, an echelle would be much more sensitive. 

Instrument builders might want to consider the advantages of CAFE in future 
spectrograph designs. From our observations alone we cannot say much about the stability 
of such conventional radial velocities over periods of years but it will be interesting to follow 
it up. 

Conventional radial velocities even of the precision reported here are not a substitute 
for those made with captive-gas fiducials in serious Doppler planet searches but the 
precision is important for observations made in spectral regions remote from telluric lines 
or captive-gas fiducials such as we have described here. 

Research funding from the Canadian Natural Sciences and Engineering Research 
Council (G.A.H.W. & E.S.) and the National Research Council of Canada (D.A.B.) is 
gratefully acknowledged. We are also indebted to the CFHT staff for their care in setting 
up the CAFE fiber system and the Gecko spectrograph. 
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Fig. 1. — A single spectrum of r Ceti flat-fielded as described in the text. The arrows define 
the 3942-3963 A region bounded by two strong lines of Al which was used in measuring the 
differential radial velocities (ARV) discussed in the paper. 
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Fig. 2. — The differential radial velocities for the five 51 Peg stars are plotted as a function 
of phase, and the unvarying star, r Ceti, as a function of time. The sine curves have the 
published planetary orbital period and K max for each star (see Table 1) and have been shifted 
in phase and ARV to give the best fit to the ARVs. For r Ceti there is a best fitting line of 
constant velocity. Below each velocity curve are shown the residuals of the ARV from the 
curves. 
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Fig. 3. — The shifts expressed as ARV determined for the Thorium Argon comparison arc 
lines as described in the text for each night. Note that the velocity scales on nights 3 and 5 
are an order of magnitude greater than on the other three nights. The dashed lines delineate 
trends. 



